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A gas bubble trapped in water by an oscillating acous- 
tic field is expected to either shrink or grow on a difltusive 
timescale, depending on the forcing strength and the bub- 
ble size. At high ambient gas concentration this has long 
been observed in experiments. However, recent sonolumines- 
cence experiments show that in certain circumstances when 
the ambient gas concentration is low the bubble can be stable 
for days. This paper presents mechanisms leading to stability 
which predict parameter dependences in agreement with the 
sonoluminescence experiments. 



Recent experiments on sonoluminescence (SL) [Q-Q al- 
low detailed studies of the dynamics of a bubble levitated 
in a periodically modulated acoustic field. Besides the 
light emission itself, one of the greatest mysteries is how 
the bubble can exist in a stable state for many billions 
of cycles. Measurements of the time between successive 
light flashes show that the total mass of the bubble re- 
mains constant to high accuracy This result con- 
tradicts classical notions about the dynamics of periodi- 
cally forced bubbles: An unforced bubble of ambient ra- 



dius Rq dissolves over a diffusive timescale, r 



D(co- 



where po is the ambient gas density in the bubble, D 
is the diffusion constant of the gas in the liquid, cq is the 
saturated concentration of the gas in the liquid, and Coo 
is the concentration of gas in the liquid far from the bub- 
ble. A strongly forced bubble grows by rectified diffusion, 
as first discovered by Blake |§,|ll| ■ This is because when 
the bubble radius is large, the gas pressure in the bubble 
is small, resulting in a strong mass fiux into the bub- 
ble. Conversely, when the bubble radius is small there is 
a strong mass outflux. Since the diffusive time scale is 
much larger than the very short time the bubble spends 
at small radii, gas cannot escape from the bubble dur- 
ing the compression phase and will be recollected during 
expansion, so that the net effect is bubble growth. At 
a special value of the ambient radius i?g rectified diffu- 
sion and normal diffusion exactly balance. However the 
above arguments suggest that this equilibrium point is 
unstable] if the ambient radius is infinitesimally different 
from i?Q the bubble is pushed away from equilibrium. 

The classic papers on rectified diffusion (see e.g. EUer 
and Crum [nO^E2|) verified the qualitative picture de- 



scribed above when Cqo/cq ~ 1. Anomalies between the- 
ory and experiment do however exist: of special note is 
Filer's observation of a stable oscillating bubble per- 
sisting over long periods. 

There are two controlled parameters in the SL experi- 
ments: the forcing pressure Pa and the gas concentration 
Coo. The key to the discovery of stable single bubble SL 
(whose existence completely contradicts the above sce- 
nario) by Gaitan et al. was that (i) Coo/co <C 1, and (n) 
Pa must lie between a lower critical pressure w l.latm 
and an upper critical pressure « l.Satm. When Pa is 
within this window of stability, the ambient radius re- 
mains constant for billions of cycles, as evidenced by the 
constant phase (p of the light emission relative to the os- 
cillatory forcing. Barber and Putterman showed that 
the "jiggle" in the phase differs by less than 50 picosec- 
onds from cycle to cycle. Outside the window of stabil- 
ity, (j) (and hence i?o) varies on a diffusive timescale [|]. 
The phase grows (implying growth of the ambient ra- 
dius), until the bubble becomes parametrically unstable 
p3[ and microbubbles pinch off. Fxperiments Q show 
that (p can oscillate indefinitely on a diffusive timescale 
via this mechanism (diffusive growth followed by pinch- 
ing off a microbubble). 

The dependence of stable SL on the gas concentration 
Coo is demonstrated by the UCLA experiments on pure 
argon bubbles For Coo/cq between approximately 

0.06 and 0.25, (f) oscillates on a diffusive timescale as de- 
scribed above. At lower argon concentration Coo/co = 
0.004 however the phase becomes perfectly stable 0, in- 
dicating a stable equilibrium. 

The goal of this paper is to suggest mechanisms lead- 
ing to stabilization: When Coo/co is decreased at forcing 
pressures Pa > l.latm, the classical unstable equilbrium 
point i?Q undergoes an inverse pitchfork bifurcations and 
actually stabilizes. At even higher forcing pressures, there 
can be several stable fixed points, although far from the 
equilibrium point small bubbles shrink and large bub- 
bles grow. This mechanism is sufficient to explain the 
stable equilibria in the SL experiments. However, high 
pressures and temperatures within the bubble cause non- 
diffusive effects [Q which can also stabilize the bubble. 
Both mechanisms predict parameter dependences con- 
sistent with SL experiments. We suggest that the dis- 
cretization of the ambient radius predicted by the diffu- 
sive mechanism provides a clear experimental signature 
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as to which effect is primarily responsible for stable SL. 

We first set up a formalism for studying the stability 
of the equilibrium point, following Fyrillas and Szeri Q 
and Lofstedt et al. 0. Let c{r,t) denote the concen- 
tration of gas dissolved in the liquid a distance r from 
the center of the bubble. For r > R{t), where R{t) is 
the radius of the bubble, c satisfies a convection diffusion 
equation 



R^R 



drC ^ DV'^C. 



(1) 



The boundary conditions are given by Henry's law 
c{R,t) — CQP{R,t)/PQ andbyc(oo,t) = Coo- The concen- 
tration gradient at the boundary gives the mass loss/gain 
of the bubble M = ATTR^Ddrc\R(^t)- 

These equations determine the growth of the bubble as 
a function of time. There are two crucial observations: 
First, EUcr noted that changing coordinates ioh — ^ ~^ 

and T = /* R'^dt transfers equation (|^) to the simpler 
form 



drC = Ddh 1 



R^ 



4/3 



dhc = 0. 



(2) 



For the following it is convenient to define the r-average 
of a function f{t) by {f{t))r = J f{t)R{t)^dt/ J R{t)^dt. 

The second observation |lJJ^ is that the bubble radius 
changes over a much faster timescale than the ambient ra- 
dius. By averaging equation ^ over the fast time scale, 
dhc{T) can be computed on the slower diffusive timescale. 
Then the dynamics of the ambient radius is given by 



R; 



dRo 



= D 



{C)r 



Jo 



dh 



(3) 
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Equilibrium points satisfy 

{P)r 



Po 



Co 



(4) 
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The equilibrium is stable if the quantity /3 — 
positive. 

Now we proceed to analyze this model. We calculate 
numerically {p)r as a function of i?o, for different driving 
pressures, by r averaging solutions R{t) of the Rayleigh 
- Plesset (RP) equation. The RP equation [^15 1^ gov- 
erns the dynamics of an acoustically forced bubble, and 
is given by 



RR + 



-(p(i?,t)-P(t)-Po) 

R d , . ^^ . R 2ct 

{p{R,t)~P{t))~^v- 



PwCw dt R PwR 

We use parameters corresponding to ||r^,||,^ an air bub- 
ble in water: the surface tension of the air-water in- 
terface is (J = O-OTSkg/s"^, while water has viscosity 



(5) 



ly = 10~®TO^/s, density p^, = lOOOfcg/m^, and speed 
of sound = 1481m/s. The acoustic field is driven 
via P{t) — PaCos{LLjt) with uj/2tt — 26AkHz and exter- 
nal pressure Pq = latm. The pressure inside the bub- 
ble varies adiabatically hke p{R) ^ (i?^ - a^)-! -*. Here 
a = i?o/8.73 is the hard core van der Waals radius. 

Figure 1 shows {p)t/Po for several values of Pa- For 
small Pa , {p)t monotonically decays with Rq, signaling 
a diffusively unstable equilibrium. For example, when 
Coo/co — 1 with a forcing amplitude of Pa = 0.8 the 
unstable equilibrium occurs at Rq « 5p,m. Note that for 
large Rq the bubble may become unstable with respect 
to shape oscillations ||l3|] . 

At large Pa however {p)t- develops oscillations as a 
function of Rq, so for a range of Coo/cq there are sev- 
eral stable equilibrium points. These stable equilibria 
only occur at low Cqo/cq, immediately suggesting a rea- 
son why diffusively stable SL only occurs under these 
conditions. As an example, see the inset of figure 1: 
when Coo/co = 10^^ and Pa = 1.25atm, there are sta- 
ble equilibria (denoted by small dots in the figure) at 
Ro = 6.5,6.8,7.1,7.5,8.0 and 8.5/im. To further ver- 
ify the existence of multiple stable equilibria, we have 
solved the full equations (|^) and (|h numerically with a 
standard finite difference scheme [ll8|. Figure ^ shows 
the ambient radius as a function of time for two different 
initial conditions with Pa = 1.25atm. In each case, the 
ambient radius saturates towards a constant (7.1/im and 
8.4/im) at long times. 

We now outline in detail the predictions of these calcu- 
lations for the SL experiments. A standard experimental 
protocol is to slowly increase the driving pressure 
Pa- The initial ambient radius depends on the prepa- 
ration of the bubble. At low pressures, if the bubble 
size is below the diffusive equilibrium curve sketched in 
figure [?], the bubble shrinks. As the forcing pressure 
is increased, there is a critical pressure where the bub- 
ble size becomes greater than Rq; the calculations for 
Coo/cq — 0.25 indicate this occurs near latm, in accord 
with experiments Above this forcing pressure, the 
bubble grows by rectified diffusion. When the ambient 
radius becomes too large the bubble is parametrically un- 
stable; in experiments the bubble decreases its radius by 
pinching off a microbubble (indicated by the downward 
arrows in figure 2). As the forcing pressure is further 
increased, the bubble size tracks the parametric insta- 
bility line |l3j. However, at the pressure Pa ~ l.latm 
stable equilibrium points appear. A sketch of this in- 
verted pitchfork bifurcation in the Rq — Pa phase dia- 
gram is shown in figure |^. If the bubble is attracted to 
one of the stable points the radius stabilizes, causing a 
discontinuous jump in the ambient radius, as observed by 
Barber et. al. Q for air bubbles. This analysis predicts 
a similar jump also for argon bubbles. The stable state 
persists for 1.1 < Pa < 1.3; Above Pa = 1.3 the equi- 
librium point destabilizes, and the bubble must return 
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to diffusive growth followed by microbubble pinching to 
survive. For even larger Pa the bubble becomes unsta- 
ble with respect to shape oscillations This entire 
scenario suggested by the classical equations of diffusive 
dynamics is in good agreement with the findings from the 
SL experiments. 

The stable equilibrium points are ultimately due to 
oscillations in (p)t as a function of Rq, which arise from 
resonances in the Rayleigh-Plesset equation. Oscillations 
even occur in the maximum radius as a function of Rq, 
so that in some situations adding more gas to a bub- 
ble decreases its maximum size. A comparison with the 
Mathieu equation is instructive: If the eigenfrequency (in 
the RP equation this depends on Rq) is an integer or half 
integer fraction of the forcing frequency, the amplitude 
of the oscillations is anomolously large. 

Our predictions are in qualitative agreement with the 
experiment; precise quantitative agreement requires ac- 
counting for several neglected effects. These include real- 
istic heat transfer [p^|6| and equations of state for the gas 
pO| , as well as shocks within the bubble |21 2^, which 
cause the light emission itself. 

Perhaps more importantly, the high pressures and tem- 
peratures within the bubble cause complications which 
influence the position and stability of the equilibria. For 
example, the RP equation predicts that the minimum 
pressure inside the bubble can become as low as 10^'^atm. 
However, the presure cannot fall below the equilibrium 
vapor pressure of water around 10~^atm. Thus, the RP 
equation can drastically underestimate the pressure in- 
side the bubble near its maximum radius. Calculations 
where the pressure inside the bubble is given by rnaa;(0.05 
atm,P(i?(t))) show that this shifts the position of the 
equilibrium point to larger Rq , without affecting the sta- 
bility. 

The increase of the mass diffusion constant with in- 
creasing temperature and pressure affects both the stabil- 
ity and the location of the equilibrium point: When the 
pressure and the temperature inside the bubble is high, 
the diffusion constant near the bubble wall is larger than 
the diffusion constant in the bulk liquid. This results in 
nondiffusive mass ejection when the bubble is small 0. 

We study this increase in the interfacial diffusion con- 
stant with an extremely simple model: Whenever the 
pressure inside the bubble exceeds a critical pressure 
Pthres{R) we discoutinuously increase the diffusion con- 
stant near the bubble wall by a factor fthres- The diffu- 
sion constant in the bulk liquid remains constant, since 
high pressures and temperatures are localized near the 
bubble wall. 

Numerical simulations of the full equations with 
and without this effect demonstrate that the unstable 
equilibrium point can be stabilized by this effect. The 
position of the equilibrium point is in general shifted to 
larger radii. Although it is difficult to determine the pre- 
cise parameter dependence of this mechanism without a 



more accurate model for the diffusion constant, it clearly 
only operates at the high pressures where sonolumines- 
cence occurs. 

To summarize, we have presented two different mech- 
anisms leading to stabilization of the ambient bubble ra- 
dius. Contrary to classical intuition, the diffusive dynam- 
ics themselves have multiple stable equilibrium points in 
the SL regime. The dependence of the mass diffusion 
constant on temperature and pressure also leads to sta- 
bilization. Which effect dominates the SL experiments? 
We suggest that this issue can be settled experimentally 
by determining whether the ambient radius takes on only 
a discrete set of values, as predicted by the purely diffu- 
sive mechanism. 
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FIG. 1. {p)t/Po as a function of Ro (in fim) 

for Pa = 0.8atm, 0.9atm,l.atm,l.latm, 1.2atm and 
1.25atm, top to bottom. Equilibrium corresponds to 
{p}t/Po = Coo/co- The equilibrium is diffusively stable if 
the slope /? = dipjr/dRo is positive. Inset: An enlarge- 
ment of Pa — 1.25atm. The straight line corresponds to 
Coo/co ~ 10^^. The intersection of the straight line with 
the curve correspond to equilibrium points. When /3 > (the 
solid dots in the figure) the equilibrium is stable. 

FIG. 2. The ambient radius Ro{t) as a function of time, 
for Pa = 1.25atm and Coo/co = 0.01, with two different ini- 
tial ambient radii. Each bubble approaches a different equi- 
librium ambient radius, demonstrating that more than one 
stable equilibrium may exist. 

FIG. 3. Sketch of the stability diagram in the _Ro — Pa 
phase space. Upon increasing the forcing Pa, a window of 
diffusive stability develops through an inverse pitchfork bi- 
furcation. Further bifurcations can occur. At larger Pa the 
stability window closes. Stable and unstable branches are 
marked by s and u. The upper curve shows the parametric 
instability line. The long arrows sketch growth by rectified 
diffusion followed by micro bubble pinchoff. 
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